Prioritizing COVID-19 vaccination efforts and dose allocation within Madagascar (preprint)

BackgroundWhile mass COVID-19 vaccination programs are underway in high-income countries, limited availability of doses has resulted in few vaccines administered in low and middle income countries (LMICs). The COVID-19 Vaccines Global Access (COVAX) is a WHO-led initiative to promote vaccine access equity to LMICs and is providing many of the doses available in these settings. However, initial doses are limited and countries, such as Madagascar, need to develop prioritization schemes to maximize the benefits of vaccination with very limited supplies. There is some consensus that dose deployment should initially target health care workers, and those who are more vulnerable including older individuals. However, questions of geographic deployment remain, in particular associated with limits around vaccine access and delivery capacity in underserved communities, for example in rural areas that may also include substantial proportions of the population. MethodsTo address these questions, we developed a mathematical model of SARS-CoV-2 transmission dynamics and simulated various vaccination allocation strategies for Madagascar. Simulated strategies were based on a number of possible geographical prioritization schemes, testing sensitivity to initial susceptibility in the population, and evaluating the potential of tests for previous infection. ResultsUsing cumulative deaths due to COVID-19 as the main outcome of interest, our results indicate that distributing the number of vaccine doses according to the number of elderly living in the region or according to the population size results in a greater reduction of mortality compared to distributing doses based on the reported number of cases and deaths. The benefits of vaccination strategies are diminished if the burden (and thus accumulated immunity) has been greatest in the most populous regions, but the overall strategy ranking remains comparable. If rapid tests for prior immunity may be swiftly and effectively delivered, there is potential for considerable gain in mortality averted, but considering delivery limitations modulates this. ConclusionAt a subnational scale, our results support the strategy adopted by the COVAX initiative at a global scale.

Existing human mobility data sources poorly predicted the spatial spread of SARS-CoV-2 in Madagascar (preprint)

For emerging epidemics such as the COVID-19 pandemic, quantifying travel is a key component of developing accurate predictive models of disease spread to inform public health planning. However, in many LMICs, traditional data sets on travel such as commuting surveys as well as non-traditional sources such as mobile phone data are lacking, or, where available, have only rarely been leveraged by the public health community. Evaluating the accuracy of available data to measure transmission-relevant travel may be further hampered by limited reporting of suspected and laboratory confirmed infections. Here, we leverage case data collected as part of a COVID-19 dashboard collated via daily reports from the Malagasy authorities on reported cases of SARS-CoV-2 across the 22 regions of Madagascar. We compare the order of the timing of when cases were reported with predictions from a SARS-CoV-2 metapopulation model of Madagascar informed using various measures of connectivity including a gravity model based on different measures of distance, Internal Migration Flow data, and mobile phone data. Overall, the models based on mobile phone connectivity and the gravity-based on Euclidean distance best predicted the observed spread. The ranks of the regions most remote from the capital were more difficult to predict but interestingly, regions where the mobile phone connectivity model was more accurate differed from those where the gravity model was most accurate. This suggests that there may be additional features of mobility or connectivity that were consistently underestimated using all approaches, but are epidemiologically relevant. This work highlights the importance of data availability and strengthening collaboration among different institutions with access to critical data - models are only as good as the data that they use, so building towards effective data-sharing pipelines is essential.

Cross-sectional cycle threshold values reflect epidemic dynamics of COVID-19 in Madagascar (preprint)

As the national reference laboratory for febrile illness in Madagascar, we processed samples from the first epidemic wave of COVID-19, between March and September 2020. We fit generalized additive models to cycle threshold (Ct) value data from our RT-qPCR platform, demonstrating a peak in high viral load, low-Ct value infections temporally coincident with peak epidemic growth rates estimated in real time from publicly-reported incidence data and retrospectively from our own laboratory testing data across three administrative regions. We additionally demonstrate a statistically significant effect of duration of time since infection onset on Ct value, suggesting that Ct value can be used as a biomarker of the stage at which an individual is sampled in the course of an infection trajectory. As an extension, the population level Ct distribution at a given timepoint can be used to estimate population-level epidemiological dynamics. We illustrate this concept by adopting a recently-developed, nested modeling approach, embedding a within-host viral kinetics model within a population-level Susceptible-Exposed-Infectious-Recovered (SEIR) framework, to mechanistically estimate epidemic growth rates from cross-sectional Ct distributions across three regions in Madagascar. We find that Ct-derived epidemic growth estimates slightly precede those derived from incidence data across the first epidemic wave, suggesting delays in surveillance and case reporting. Our findings indicate that public reporting of Ct values could offer an important resource for epidemiological inference in low surveillance settings, enabling forecasts of impending incidence peaks in regions with limited case reporting.

Monitoring for outbreak associated excess mortality in an African city: Detection limits in Antananarivo, Madagascar (preprint)

Quantitative estimates of the impact of infectious disease outbreaks are required to develop measured policy responses. In many low- and middle-income countries, inadequate surveillance and incompleteness of death registration are important barriers. Here, we characterize how large an impact on mortality would have to be to be detectable using the uniquely detailed mortality notification data from the city of Antananarivo in Madagascar, with application to a recent measles outbreak. The weekly mortality rate of children during the 2018-2019 measles outbreak was 154% above the expected value at its peak, and the signal can be detected earlier in children than in the general population. This approach to detecting anomalies from expected baseline mortality allows us to delineate the prevalence of COVID-19 at which excess mortality would be detectable with the existing death notification system in the capital of Madagascar. Given current age-specific estimates of the COVID-19 fatality ratio and the age structure of the population in Antananarivo, we estimate that as few as 11 deaths per week in the 60-70 years age group (corresponding to an infection rate of approximately 1%) would detectably exceed the baseline. Data from 2020 will undergo necessary processing and quality control in the coming months. Our results provide a baseline for interpreting this information.

High variation expected in the pace and burden of SARS-CoV-2 outbreaks across sub-Saharan Africa (preprint)

A surprising feature of the SARS-CoV-2 pandemic to date is the low burdens reported in sub-Saharan Africa (SSA) countries relative to other global regions. Potential explanations (e.g., warmer environments [1], younger populations [2,3,4]) have yet to be framed within a comprehensive analysis accounting for factors that may offset the effects of climate and demography. Here, we synthesize factors hypothesized to shape the pace of this pandemic and its burden as it moves across SSA, encompassing demographic, comorbidity, climatic, healthcare and intervention capacity, and human mobility dimensions of risk. We find large scale diversity in probable drivers, such that outcomes are likely to be highly variable among SSA countries. While simulation shows that extensive climatic variation among SSA population centers has little effect on early outbreak trajectories, heterogeneity in connectivity is likely to play a large role in shaping the pace of viral spread. The prolonged, asynchronous outbreaks expected in weakly connected settings may result in extended stress to health systems. In addition, the observed variability in comorbidities and access to care will likely modulate the severity of infection: We show that even small shifts in the infection fatality ratio towards younger ages, which are likely in high risk settings, can eliminate the protective effect of younger populations. We highlight countries with elevated risk of slow pace, high burden outbreaks. Empirical data on the spatial extent of outbreaks within SSA countries, their patterns in severity over age, and the relationship between epidemic pace and health system disruptions are urgently needed to guide efforts to mitigate the high burden scenarios explored here.